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Abstract 

Relational particle mechanics is useful for modelling whole-universe issues such as quantum cosmology or the problem 
of time in quantum gravity, including some aspects outside the reach of comparably complex minisuperspace models. 
In this article, we consider the mechanics of pure shape and not scale of 4 particles on a line, so that the only physically 
significant quantities are ratios of relative separations between the constituents' physical objects. Many of our ideas 
and workings extend to the N-particle case. As such models' configurations resemble depictions of metro lines in public 
transport maps, we term them 'N-stop metrolands'. This 4-stop model's configuration space is a 2-sphere, from which 
our metroland mechanics interpretation is via the 'cubic' tessellation. This model yields conserved quantities which are 
mathematically SO(3) objects like angular momenta but are physically relative dilational momenta (i.e. coordinates 
dotted with momenta). We provide and interpret various exact and approximate classical and quantum solutions 
for 4-stop metroland; from these results one can construct expectations and spreads of shape operators that admit 
interpretations as relative sizes and the 'homogeneity of the model universe's contents', and also objects of significance 
for the problem of time in quantum gravity (e.g. in the naive Schrodinger and records theory timeless approaches) . 
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1 Introduction 



Euclidean relational particle mechanics (ERPM) (proposed in [J and further studied in pil51Hl[51[51[71igi HUl[TTin2irn] 
is a mechanics in which only relative times, relative angles and relative separations are meaningful. On the other hand, in 
similarity relational particle mechanics (SRPM) (proposed in [TS] and further studied in [T5J |SJ [T71 171 [T5J |HJ [1511201 US]), only 
relative times, relative angles and ratios of relative separations are meaningful. More precisely, these theories implement 
the following two Barbour- type relational postulates. 

1) They are temporally relational [U 121] 1221 123} 118] . i.e. there is no meaningful primary notion of time for the whole 
system thereby described (e.g. the universe), which is implemented by using actions that are manifestly reparametrization 
invariant while also being free of extraneous time-related variables [such as Newtonian time or General Relativity (GR)'s 
lapse] . This reparametrization invariance then directly produces primary constraints quadratic in the momenta. 

2) They are configurationally relational, which can be conceived in terms of a certain group G of transformations that 
act on the theory's configuration space Q being held to be physically meaningless [U [211 E21 [23l HH [13]. This can be 
implemented by such as using arbitrary-G-frame-corrected quantities rather than 'bare' Q-configurations. For, despite this 
augmenting Q to the principal bundle P(Q,G), variation with respect to each adjoined independent auxiliary G- variable 
produces a secondary constraint linear in the momenta which removes one G degree of freedom and one redundant degree 
of freedom among the Q variables. Thus, one ends up dealing with the desired reduced configuration space - the quotient 
space Q/G. Configurational rclationalism includes as subcases both spatial rclationalism (for spatial transformations) 
and internal relationalism (in the sense of gauge theory) . 

For ERPM, the Jacobi-type (Mj action isQ 

I =2 J d\y/T{E-V} , with T = i - k x R l } 2 /2 , (1) 

and for SRPM, our presentation of it is 

I = 2 /" d\y/T{E-V} , withT = ^\ x R l + cR i } 2 /2I . (2) 

These implement the above relational postulates for the corresponding Euclidean and similarity G's [TJ [15l [6l [TTJ [9]. 
Equivalent theories formulated directly in terms of rotational (and dilational) invariant quantities can also be arrived at 
by considering the space of shapes and mechanics thereupon |30l 118] . It is then of interest what structure one gets when 
one quantizes such theories El [HI [111 [M 123 1101 [121 [131 [111 I2H] • 

The Barbour-type indirect formulation of RPM's ( |1|2[ ), moreover, is particularly interesting through how the geometro- 
dynamical form of GR can be cast in direct parallel: it also obeys postulates 1) and 2) implemented as follows [2"T1 12"3"1 l3"T] 
(some features of which are already anticipated in [32] ) 

S GR = 2 J d\ J d 3 xVh^T GR {mc(h) - 2A} for T GR = ^M» vpa - £fih iw }{h pa - £ph p<J } ; (3) 

in this case, Q is the space Riem(E) of Riemannian 3-metrics on a fixed spatial topology S, and G is the corresponding 
3-diffeomorphism group, Diff(E). 



The way that the physical equations follow from each of the above actions then has many parallels. By reparametriza- 
tion invariance |34] each has a primary constraint quadratic in the momenta: for GR the Hamiltonian constraint 

U = A/^tt" V - y/h{Bic(h) - 2A} = (4) 

'RPM's are relational in Barbour's sense of the word rather than Rovelli's distinct one; see e.g. |591[3l[4] for these authors' original material 
and [9] for a discussion of some differences. 

2 R l are relative Jacobi coordinates 1251 : linear combinations of relative particle separation vectors that produce a diagonal kinetic term and 
are particular inter-particle cluster separation vectors with associated cluster masses Lower-case Latin indices run over 1 to n = N — 1 for 
N the number of particles, and lower-case greek ones are spatial indices; the spatial dimension is d. A is label time and ' is the derivative with 
respect to this. Using such relative coordinates, one has already incorporated the highly trivial translation part of the Euclidean or similarity 
groups, b is a rotational auxiliary velocity whereby the rotation part of these groups is implemented. In the SRPM case, c is a dilational 
auxiliary velocity implementing the additional scaling part. Mi a jp = ln&ijS a p is the mass matrix with determinant M and inverse N la ^ . P t 
is the momentum conjugate to W . I is the moment of inertia, /J,i\R l \ 2 . For ERPM, T, V and E are kinetic, potential and total energy 
terms with the usual physical dimensions. In our 'pure shape' formulation of the SRPM, the kinetic term T has dimensions of (energy)// and 
E — V has dimensions of (energy) X 7. Consistency dictates that this V additionally be a homogeneous function of the W; in fact, in the given 
'pure shape' formulation, it must be homogeneous of degree zero. This turns out not to be a heavy restriction due to I being constant after 
variation and useable to homogenize (see Sec 2 for an example) . While an actual energy is prohibited by this consistency, the abovementioncd 
constant E is permissible instead. 

3 The spatial topology S is taken to be compact without boundary. h pv is a spatial 3-metric thereupon, with determinant h, covariant 
derivative D^, Ricci scalar Ric(/i) and conjugate momentum ir^" . A is the cosmological constant. J\4^" pa = h^ p h v<T — h pu h pa is the inverse 
DeWitt supermetric with determinant A4 and inverse N^vpo- To represent this as a configuration space metric (i.e. with just two indices, and 
downstairs), use DeWitt's 2 index to 1 index map 1331 . is the velocity of the frame; in the manifestly relational formulation of GR, this 
cyclic velocity plays the role more usually played by the shift Lagrange multiplier coordinate. £p is the Lie derivative with respect to . 
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and, for ERPM and SRPM respectively, the 'energy constraints': 

q = N ia ^P ia P j0 /2 + V = E , q = IN ia ^P ia P j0 /2 + V = E . (5) 

By variation with respect to the auxiliary G- variables, each relational theory has constraints linear in the momenta: for 
GR, the momentum constraint 

= -2D v ir\ = (6) 

from variation with respect to F^, and, for RPM's, the zero total angular momentum and zero total dilational momentum 
constraints 

^E i 5 i x£ j =0 , D=£ j ^. J P j = (7) 

from variation with respect to & M and c (so the latter constraint occurs only in SRPM). The zero total dilational momentum 
constraint moreover closely parallels the well-known GR maximal slicing condition |35j . h^Tr^ = 0. Furthermore, 
much like generalizing maximal slicing to constant mean curvature slicing |36j turns on a 'York time' variable |37l 138] 
t York = 2h^ u /Zy r h = c(A alone), one can think of the passage from SRPM to ERPM as involving an extra 'Euler time' 
variable £ Eulcr oc ■ P4 — C(A). This is all underlied for both GR and RPM's by shape-scale splits, the role of scale 

being played by \fl or / for RPM's and by such as the scalefactor a or \fh in GR. In both cases it is then tempting to 
use the singled-out scale as a time variable but this runs into monotonicity problems which are avoided by using as times 
the quantities conjugate to (a function of) the scale, i.e. t York and t Euler . 

There are further analogies at the configuration space level. If 1) R(N,d) the relative space of relative interparticle 
(cluster) separation vectors and Riem(S) are held to be analogous, then so are 2) Relational space = R(N, c?)/Rot(c?) 
for Rot(d) the d-dimensional rotations and superspace(E) = Ricm(E)/Diff(£). 3) Shape space — R(N, d)/Rot x Dil for 
Dil the dilations and conformal superspace [351110] CS(S) = Riem(£)/DifT(£) x Conf(E) for Conf(E) the conformal 
transformations on S. 4) The cone representation of relational space in shape-scale split variables [TT] and 
{CS + V}(£) = Riem(E)/Diff(E) x VPConf(E) [3J] for VPConf(E) the conformal transformations that preserve the 
volume of the universe, V [4"T| . Also, both these GR and RPM configuration spaces arc in general stratified, and both 
have physically significant bad points (e.g. a = is the Big Bang and / = is the maximal collision). 

There are yet more analogies [HI [331 IHZ1 H31 IIH EH S US HH1 HI 113 H7J HH] at the level of various strategies toward the 
resolution of the Problem of Timc|^] and various other aspects of quantum cosmology. The above quadratic constraints 
give frozen (i.e. timeless or stationary) quantum equations. For GR, this is the Wheeler-DeWitt equation, 

= -ft 2 ' {^^7 {V^A^I^} - eRic(M)*} ' - Vh(mc(h) - 2A}* = , (8) 

where "J is the wavefunction of the universe; ' ' implies in general various well-definedness issues (see e.g. |19j for a sum- 
mary) and need for a choice of operator-ordering (we use conformal ordering in this paper, c.f. Sec 3.1). Correspondingly, 
for RPM's, 

^ = l-^^K {n ab Vm^X -£Ric(M)}y + \/y = EV . (9) 



For the moment a = ia, B = j(3 and the Q A are the R la . A AB is the inverse mass matrix for ERPM and / times it for 
SRPM. However, we use this equation more generally than that below for reduced RPM's in which linear constraints have 
been taken care of, this being explicitly possible in 1- or 2-d [T5] , 

An important feature of GR (and one missed out by minisuperspace models [49, 50] ) is that of linear constraints causing 
substantial complications e.g. in attempted resolutions of the Problem of Time. That is the momentum constraint for GR 
or its quantum counterpart], while RPM's have the linear constraints or their quantum counterparts]. However, 
minisupcrspaces, unlike RPM's, have more specific and GR-inhcritcd potentials and indefinite kinetic terms. Thus both 
minisuperspace and RPM's are valuable in complementary ways as toy models^] 

Some of the strategies toward resolving the Problem of Time are as follows. 

A) Perhaps one is to find a time hidden within classical GR |38j and thus obtain a wave equation that depends on it from 
the outset at the quantum level. York time is a GR example of such and Euler time is an ERPM model of it. 

B) Perhaps one has slow, heavy 'H' variables that provide an approximate timestandard with respect to which the 
other fast, light 'L' degrees of freedom evolve [52J [38], [47] . In quantum cosmology the role of H is played by scale (and 
homogeneous matter modes), so ERPM's in scale-shape split are more faithful scmiclassical models of this than SRPM's 
themselves can muster. 

C) A number of approaches take timelessness at face value. One considers only questions about the universe 'being', 
rather than 'becoming', a certain way. This can cause at least some practical limitations, but nevertheless can address 



4 This notorious problem occurs because 'time' takes a different meaning in each of GR and ordinary quantum theory. This incompatibility 
underscores a number of problems with trying to replace these two branches with a single framework in situations in which the premises of 
both apply, namely in black holes and in the very early universe. One facet of the Problem of Time that shows up in attempting canonical 
quantization is that the lack of linear momentum dependence of the GR Hamiltonian constraint, but it does have many other facets |38| . 

5 Midisuperspace 1511 unites all these desirable features but is unfortunately then calculationally too hard for many of the strategies. 
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at least some questions of interest. For example, Hawking and Page's naive Schrddinger interpretation [53] concerns the 
'being' probabilities for universe properties such as: what is the probability that the universe is large? Flat? Isotropic? 
Homogeneous? One obtains these via consideration of integrals of |W| 2 over suitable regions of the configuration space. 
This approach is termed 'naive' due to it not using any further features of the constraint equations. The conditional 
probabilities interpretation [43 goes further by addressing conditioned questions of 'being' such as 'what is the probability 
that the universe is flat given that it is isotropic'? Records theory [43j [55j |45l SJ [56j [57] involves localized subconfigurations 
of a single instant - whether these contain useable information, are correlated to each other, and whether a semblance 
of dynamics or history arises from this. RPM's are superior to minisuperspace for such a study as they have a notion 
of localization in space, and more options for well-characterized localization in configuration space through their kinetic 
terms possessing positive-definite metrics. 

D) Perhaps instead it is the histories that are primary {histories theory [55 , 58 ). There is a records theory within histories 
theory, and histories decohereing is one possible way of obtaining a semiclassical regime in the first place, making B) to 

D) of particular interest to one of us [571 ITU] . 

E) Distinct timeless approaches involve evolving constants of the motion ('Heisenberg' rather than 'Schrodinger' style 
QM), or partial observables |59j (which are used in Loop Quantum Gravity's master constraint program [60 ). 

Some approaches to the Problem of Time that do not have an RPM analogue include superspace time (which requires 
indefinite configuration spaces) and third quantization (which requires field theoretic rather than finite models). 

We are in the process of building up a reasonable set of RPM models, paralleling e.g. the development of minisuperspace 
in the early 70's [H0[5U], or Carlip's work in the 90's for 2+1 gravity [HT] (see [35] for yet further useful toy model arenas 
for Problem of Time approaches). Also RPM's serve as a bridge from highly-studied ideas in molecular physics and 'mini- 
and midi'superspace, which may serve to import ideas and tools from the former to the latter. 

Our build-up is for RPM's in 1-d and 2-d; for N particles, we term these, respectively, N-stop metroland and N-a- 
gonland (the first two nontrivial N-a-gonlands we furthermore refer to as triangleland and quadrilateralland) . We choose 
to study these models because their configuration spaces are highly tractable mathematically j3Ul[IB] : § N ~ 2 spheres in 
1-d and CP N_2 complex projective spaces in 2-d. This is for the choice of plain shapes rather than oriented shapes. 
[I.e. we make the choice of treating each shape and its mirror image as distinct; in this paper's 4-stop metroland model, 
this means that we regard the 1,2,3,4 ordering of the particles to be distinct from the 4,3,2,1 one. The opposite choice 
gives configuration spaces S fe /Z 2 = RP k (real projective spaces), and CP fc /Z 2 which are somewhat harder to model. 1 
N.B. that the interesting theoretical parallels between GR and RPM's are unaffected by our choice of plain shapes and 
of low-d RPM's. We are presently studying scalefree models as these are more straightforward than models with scale 
(though we will need to move on to scaled models as regards reasonably quantum-cosmologically realistic modelling of 
the semiclassical approach; note also that scalefree problems occur as a subproblems in models with scale [63] [TBI E], so 
studying these first also makes sense even from this semiclassical quantum cosmological perspective). 

This paper considers scalefree 4-stop metroland (the smallest scalefree metroland to have the nontrivialities associ- 
ated with having 2 physical degrees of freedom, so that one physical quantity can be expressed in terms of another, a 
feature necessary for records theory's correlations, while semiclassical approaches need at least one H and at least one L, 
decoherence only makes sense if one thing decoheres another, and so on) . Additionally, there are indications that 4-stop 
metroland is simpler than triangleland [TT] (which also has two physical degrees of freedom and spherical shape space by 
CP 1 = § 2 ), particularly in the cases with scale and at the quantum level, so that the present paper is useful toward how 
to subsequently deal with these other more complicated cases. Also, many of the present paper's workings readily extend 
to N-stop metroland. 4-stop metroland and triangleland are both useful preliminaries for studying quadrilateralland, 
which is the simplest RPM to exhibit a number of geometrical nontrivialities, including some relevant to Problem of Time 
approaches and some that are archetypal of 2-d problems in ways that triangleland is not. Moreover, 4-stop metroland 
itself is already suitable for study of various timeless approaches to the Problem of Time. 

In Sec 2 we begin with a classical treatment of 4-stop metroland in its reduced form. We give a tessellation of the shape 
sphere corresponding to 4-stop metroland's physical interpretation and provide useful and meaningful shape quantities for 
our study. We then study 4-stop metroland's equations of motion and its conserved quantities, among which some have 
angular momentum-like mathematics but are physically dilational rather than angular momenta. We interpret simple 
subcases of multiple harmonic oscillator- like potentials' solutions using our tessellation and shape quantities, concentrating 
on the case of two localized and well-separated subsystems that is motivated by our interest in timeless approaches. 

In Sec 3 we consider time-independent Schrodinger equations for these problems, by firstly interpreting their exact 
and asymptotic solutions against our tessellation. Secondly, we compute expectations and spreads of our shape quantities 
promoted to quantum-mechanical shape operators. Thirdly, we consider perturbations about the simplest case in which 
the 'springs' balance each other out to produce a constant potential and hence spherical harmonics mathematics. Fourthly, 
we note and apply a number of analogies with molecular physics to our study. We conclude in Sec 4, including discussion 
of how our model and slightly larger versions thereof can be used as an arena for investigation of a number of Problem 
of Time strategies - of which we provide naive Schrodinger interpretation examples - and by commenting on 'mini- and 
midi'superspace counterparts of this paper's shape operators. 
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2 4-stop metroland at the classical level 



2.1 Passage to reduced form of 4-stop metroland and useful coordinatizations of it 

The unreduced action is given by the SRPM case within eq. (2) further restricted to being in 1-d (so there are no rotations) 
and for 4 particles, and so 3 relative separations and thus 3 relative Jacobi coordinates, R l : 

1 = 2 J dA^/T{E — V} with T = ^2 3 .^ 1 fi i {R i + cR i } 2 /2I . (10) 

We will find it more convenient to deal with the subsequent physics in terms of l 1 = ^/JmR 1 the mass-weighted relative 
Jacobi coordinates [which are physically the square roots of the partial moments of inertia P — fiiR 12 (no sum)], and, 
after variation, their 'normalized counterparts n l — i % ji and N l = P/I for / the total moment of inertia and I = \fl. It 
will often be convenient to use n x , n y , n z for the components of n l . We take these Jacobi coordinates to be, in terms of 
particle position coordinates, Jacobi H-coordinates rather than Jacobi K-coordinates (Fig 1) with quantum cosmological 
and records theoretic applications in mind: two equal particle number clusters treated on the same footing, each could 
model the seed of a galaxy, or be a nontrivial record (of which we need at least 2 to consider correlations between records) . 

a) 2 3 b) 2 3 c) cluster {12} cluster{34} 

,1 ,2 



cluster 
{12} 




cluster I . L ^ l _ 

{34} 



-4= e e 



1 



l 3 



= arctan 



Fi gure 1: a) and b) explain in 1-d the origin of the names H- and K-coordinates. Using {a...c} to denote the cluster composed of particles 
a, c ordered from left to right, + is the centre of mass (COM) of cluster {12}, X is the COM of cluster {34} and T is the COM of the triple 
cluster {123}. c) What H-coordinates look like in 1-d: the H has been 'squashed'. 

Now let us perform some variational manoeuvres on the above action. It is useful to bear in mind from the outset 
that our 4-stop metroland's reduced configuration space is S 2 and we are trying to bring this out as cleanly as possible 
by removing extraneous variables and seeking for standard coordinates on this. Variation with respect to the dilational 
auxiliary c gives the dilational constraint ([7]), the Lagrangian form for which can be rearranged to 

and used to eliminate c from the action, producing (2) but with T rcc i in place of T: 

Then, via the coordinate transformation 

( ^FFTpTA 3 ) , $ = arctan (l 2 /^) , (13) 
( 12 ) becomes 

T S 2 = {6 2 + sin 2 e4 2 }/2 . (14) 

The coordinate ranges are < O < ir and < $ < 2ir, so these are geometrically the standard azimuthal and polar 
spherical angles on the unit shape space sphere S 2 . Inversely, 

i 1 = isinOcos<l> , i 2 — isinOsin<l> , i 3 — tcos0 . (15) 

Thus 4-stop metroland has t = \fl playing a (constant) radius role, and the l 1 are Cartesian coordinates in the 
'surrounding' Euclidean relational space i?(4, 1) = R 3 , subject to the on-sphere condition 

EL/' = ELM 2 = M 2 = I (constant) , or = ELiW' = 1 • ( 16 ) 

The n l are then the components of the unit Cartesian vector [(sin cos sin sin <f>, cos 0) in spherical polar coordinates] . 
This should be contrasted with the way the sphere arising in scalefree triangleland [9 being harder to deal with from the 
perspective of the 'surrounding' Euclidean relational configuration space, which is ]R 3 . Scalefree triangleland's i 1 and t 2 
are related to the Cartesian coordinates of the surrounding relational space not in the above familiar Cartesian way, but 
rather in the less straightforward 'Dragt' way [TT], corresponding to having to use not i but / as radial variable. 
The following formulae are also useful below: 

I ft TL — ft *2iTl TL 

cos© = n Z) sin6 = Jn x 2 + n y 2 , cos$ = x sin$ = y cos2$ = — =- — y — : sin2$ = — - v , 

v ^Jn x 2 + n 2 y/n x 2 + n v 2 n x 2 + n y 2 n x 2 + n. 2 

(17) 

and what is (from the geometrical perspective) a stereographic radial coordinate, 



K = tanf = y/{l-n z }/{l + n z } . (18) 



A 



2.2 Action and banal-conformal representations for this paper 



We take the Jacobi action corresponding to (14 1 



1 = 2 y"dAv/T s2 {E-V} 



(19) 



to be primary. As well as by Sec 2.1's reduction, this can be obtained by |18j considering a natural mechanics in Jacobi's 
geometrical sense [21] on the space of shapes [3U] . 

The 'banal' conformal invariance of each of the above product-type actions T — ► f2 2 T, E — V — > {E — V}/£7 2 is useful 
below in 'passing' factors between T and E — V (which we term 'picking a distinct banal representation'). The above-given 
forms of T and E — V are the geometrically-natural ones (both in the scale-invariant sense and in the sense of having 
the standard spherical metric on the shape space sphere) and mechanically-natural in the sense that E itself appears in 
them rather than E times some power of the moment of inertia. However, in some applications, (1Z, $) coordinates and 
= {1 + lZ 2 }/2 is useful; we term this the 'flat banal representation' as its T is flat, and denote it by tilde-ing. 





Figure 2: a) On the configuration space represented as a sphere, there are 8 triple collision (T) points and 6 double-double (DD) collision 
points. Each DD is attached to 4 T's, and each T to 3 T's and 3 DD's, in each case by single double collision lines. This forms a tessellation 
with 24 identical spherical isosceles triangle faces, 36 edges and 14 vertices. The T's and DD's form respectively the vertices of a cube and 
the octahaedron dual to it (dashed lines in the second and third subfigures), so that the physical interpretation has the symmetry group of 
the cube, of order 48. This is isomorphic to 54 X Z2 for S4 the permutation group of 4 objects, thus realizing the ways of labelling the 4 
particles and ascribing an orientation. See p 72-75 of 164] for mathematical discussion of this tessellation and 65 for another occurrence of it 
in mechanics. In this arrangement, the T's and DD's form 7 antipodal pairs, thus picking out 7 preferred axes. The 3 axes corresponding to 
antipodal DD pairs are related to the 3 permutations of Jacobi H-coordinates and the 4 axes corresponding to antipodal T pairs are related to 
the 4 permutations of Jacobi K-coordinates. This relation is in the sense that the poles in each case correspond to what each coordinatization 
picks out as intra-cluster coordinates both going to zero i.e. collapse of both clusters for an H or collapse of the triple cluster for a K. 
b) To make statements concerning shapes being near a DD or T - i.e. well-localized (intra-cluster distances far smaller than external distances 
to non-member particles/clusters), spherical caps © < e are useful in the corresponding spherical polar coordinate chart. In particular, with 
this paper's usual choice of axis, the polar caps are where there is both {12} and {34} localization, so following clusters {12} and {34} makes 
sense. Denote this clustering (i.e. partition into clusters) by {12,34}. The opposite notion is merged clusters for which the COM's of {12} 
and {34} are near each other so that these clusters largely overlap (which is, in a certain sense, a more 'homogeneous' universe model). This 
corresponds to belts it/2 — S < < 7r/2 + 5 around the equator. (Multi)lunes <I>o — r\ < <3? < <3>o + r) also correspond to physically meaningful 
statements. E.g. being in the bilune around the Greenwich meridian means that the {34} cluster is localized, being in the bilune around the 
'Bangladeshi' meridian perpendicular to the Greenwich one means that the {12} cluster is localized and being in the tetralune at 7r/4 to all 
of these signifies that clusters {12} and {34} are of similar size', i.e. r;-close to contents homogeneity (i.e. that the particle clusters that make 
up the model universe are, among themselves, of similar constitution). The above sort of approximate notions of shape are in the spirit of 
those used in e.g. Kendall et al. 1 3 1 . and we make use of the corresponding configuration space regions in our naive Schrodinger approach 
calculations in Sec 4.3. Moreover, note that notions of 'relative size' and 'similar contents' here in fact involve more concretely the \J mass X 
distance combination (whose squares are partial moments of inertia). 



2.3 Physical interpretation by tessellation, charts and shape operators 

The Jacobi H-coordinates in use are better-adapted for 'seeing' double-double collisions [see Fig 2a)] rather than triple 
collisions (the opposite is the case for Jacobi K-coordinates) , so that it is useful to preliminarily work out and graphically 
represent the mechanical interpretations of the various zones of our problem's configuration space of shapes: Fig 2. This 
is of considerable use below in interpreting classical trajectories (as paths upon this figure) and classical potentials and 
quantum-mechanical probability density functions (as height functions over this figure). Spherical polar coordinates about 
each axis in Fig 2 are natural for the study of the corresponding H or K structure. Thus each choice of H- or K-coordinates 
has a different natural spherical polar coordinate chart. Any two of these natural charts suffice to form an atlas for the 
sphere (each goes bad solely at its poles, where its axial angle ceases to be defined). To look extremely close to a pole, 
one can 'cartesianize' e.g. after projecting the relevant hemisphere onto the equatorial disc. 

The n l are interesting quantities with which to describe the shape of the configuration. Of course, only 2 of the 3 n l 



are independent, by the on-sphere restriction (16 1. For the example of H-coordinates that follows the {12,34} clustering 
arrangement that we follow in particular in this paper, the corresponding n z is a signed quantifier of the relative size of 
the universe from the perspective of an observer in either cluster. On the other hand, n x is a signed quantifier of the size 



of the universe from the less Copernican perspective of an observer specifically in cluster {12} [n y has the same meaning 
but for cluster {34}]. Thus we term these shape quantities RelSize(12,34), RelSize(l,2) and RelSize(3,4) respectively]^] 

RelSize(12,34) small means that clusters {12} and {34} are merged, and corresponds geometrically to the equatorial 
belt. RelSize(12,34) large means physically that clusters {12} and {34} each cover but a small portion of the model 
universe, and corresponds geometrically to the polar caps. RelSize(l,2) small means physically that cluster {12} is but 
a speck in the firmament, and corresponds geometrically to a belt around the 'Bangladeshi' meridian. RelSize(l,2) large 
means physically that cluster {12} engulfs the rest of the model universe, and corresponds geometrically to an antipodal 
pair of caps around each of the intersections of the equator and the Greenwich meridian. 

A quantifier of the contents inhomogeneity between the two clusters is which is related to the ratio of the size of 
{34} to that of {12}'s by (13 1. N.B. that the last 3 paragraphs refer, more concretely, to y / massx length, so that large 
mass hierarchies can distort intuitive notions of 'actual size'. 



2.4 Rotor and planar mechanics analogies for 4-stop metroland 

By inspection of the kinetic term, there are clear analogies between this 4-stop metroland problem and well-known rotor 
and planar problems in ordinary mechanics. For the first analogy, 

arctan(-\/RelSize(12) 2 + RclSize(34) 2 /RelSize(12,34)) = G < — > 9 = (azimuthal coordinate of the axis in space) (20) 

arctan (RelSize(34)/RelSize(12)) = <f> < — > <j> — (polar coordinate of the axis in space) , (21) 

1 *-* -frotor ( moment of inertia of the rotor ) . (22) 

For the second analogy, transform B to the radial stereographic coordinate 1Z — tan® and pass to the 'tilded' banal 
representation. One then has the flat plane polar coordinates kinetic term, so 

-RclSize(12,34)}/{l + RelSize(12, 34)} = K < — > r = (radial coordinate of test particle) , (23) 

arctan (RelSize(34)/RelSize(12)) = <f> < — > <f> — (polar coordinate of test particle) , (24) 

1 <-> m = (test particle mass) . (25) 

These analogies will be furtherly fruitful in analyzing 4-stop metroland's equations of motion and conserved quantities in 
the next 2 subsections, as well as when further specifics about the potential are brought in (see Sec 3.6). 



2.5 Equations of motion for 4-stop metroland 

The equations of motion are 

9** - sin 6 cos 6 <T 2 = -V e , {sin 2 6$*}* = -V,* . (26) 

(The star is derivative with respect to the relational approach's emergent time t: * = d/dt = a/{E- V}/T; for which 
the equations of motion simplify. This is readily deduced to banal-transform as * — > f2 -2 * |55].] 

V is independent of A itself and so one of these can be replaced by the 'energy relation' (a first integral): 

{6* 2 +sin 2 6<r 2 }/2 + V(6,$) = E , constant . (27) 

If the potential is additionally ^-independent (which we term 'special' and whose planar mechanics analogue is termed 
central), then the <1>-Euler- Lagrange equation above gives another first integral, 

sin 2 6<r = X> . (28) 

For both of the analogies above, the correpsonding SO(2) or SO(3) related constant of the motion has the physical meaning 
of an angular momentum; for its interpretation in the present context, however, see the next subsection. In the special 
case, one can now furthermore combine the last 2 equations in two ways. Firstly, 

E = 6* 2 /2 + P 2 /2sin 2 6 + V(0) = 6* 2 /2 + V cff (29) 

which in the planar central problem amounts to modification of the potential by a centrifugal barrier, while, in the more 
directly analogous rotor problem, amounts to placing a centrifugal barrier at each pole. In our problem, it takes the 
latter 'bipolar barrier' form. Secondly, (for V ^ 0) ! $ 2 /2 = sin 4 0{E — V c ff}/2? 2 . Both of these straightforwardly give 
quadratures relating to, respectively, t (orbit traversal rate) and <& (shape of the orbit, the V = case giving a <f> = 
constant l-d motion without any double barrier). 

If V is also ©-independent and thus constant, we get 3 P-quantities from freedom to pick whichever axis to have a 
conserved <& about. We call this constant-potential case the 'very special case' (the counterpart of which in the second 
analogy is the rigid rotor). 

6 There are other such operators corresponding to attaching significance to other clustering arrangements obtained by permuting the particles 
in defining the Jacobi H-coordinates, while similar quantities can be defined for the various permutations of Jacobi K-coordinatcs; see [20] for 
a brief account of the shape operators for each of these. RelSize(12,34) is also \/l — n x 2 — n H 2 , so it can also be viewed as a 'measure' of 
noncollapse of at least one of the model universe's clusters; one can readily work out such 'dual statements' for other shape quantities. 
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2.6 Further discussion of 4-stop metroland's conserved quantitites 

One interesting issue in RPM's is what these angular momentum-like quantities are physically Trianglcland is spatially 
2-d and as such affords a notion of angular momentum; its conserved quantity J turns out to be the relative angular 
momentum between its subsystems. But the present paper's 4-stop metroland problem, however, is spatially l-d, so no 
angular momentum in space (relative or otherwise) is possible. What then is the meaning of the conserved quantity T> in 
terms of the n % or RelSize variables? 

V = n x n y * - n y n x * , (30) 

which is the '3-component of an angular momentum in the Euclidean relational configuration space i?(4, 1) = M 3 '. 
Moreover, using for individual/partial dilations FtPi (no sum) 

V = D 2 n x / n y - Din y /n x (31) 

so that D is a (weighted) relative dilational quantity corresponding to a particular exchange of dilational momentum 
between the {12} and {34} clusters. 

In the very special case, there are 3 T> conserved quantities forming a vector in the 'surrounding' Euclidean relational 
configuration space ]R 3 , of which the above V is the 3-component: 

V, = e ijk n j n k * = D k n 3 /n k - D 3 n k /n 3 (32) 

where i, j, k are a cycle of 1, 2, 3. 

All in all, the less special a problem is, the more types of relative dilational momentum exchanges it has. 

That we get 'angular momentum mathematics', we explain as follows. The body of mathematics habitually associ- 
ated with angular momentum can actually be associated more generally (in terms of what physics it covers, not what 
mathematics it is, as further explained in Appendix A) with rational (i.e. 'ratio-based') quantities rather than just with 
angular ones (which are a subset thereof)]^] Therefore 'rational momentum mathematics' would appear to be a more 
widely appropriate term, covering both angular momentum and dilational momentum as subcases. The objects in ques- 
tion continue to possess antisymmetry in this more general setting since this derives from differentiating a (function of a) 
ratio by the (chain rule and) quotient rule: {f(y/x)}* — f'{xy* — yx*}/x 2 and thus occurs irrespective of whether that 
ratio admits an interpretation as an angle in physical space. 

2.7 Passage to 4-stop metroland's Hamiltonian 

The conjugate momenta arc then 

Pe = 6* , p® = sin 2 6 $*=£>. (33) 
[Also, T>i — €ij k Jpk for pk the momentum conjugate to t k .] The momenta obey a quadratic constraint 

Q = Pe 2 /2 + p$ 2 /2sin 2 e + V(9, $) = E , (34) 

the middle expression of which also serves as the classical Hamiltonian H for the system. 

2.8 Harmonic oscillator like potentials for 4-stop metroland 

With eventual timeless records and structure formation goals in mind, we intend to follow a particular clustering - the 
{12,34} one - using a particular permutation of H-coordinates which is physically picked out by considering not the 
most general array of 6 springs between the particles but rather the following. We take the mechanical picture in Jacobi 
coordinates as primary and consider springs within each of the {12} and {34} clusters and between the centres of mass 
of the two clusters (reinterpretable if one so wishes as a superposition of inter-particle springs) . Then the potential is 

V = 3 i=1 K ii ni } 2 / 2 = i K i sin20 cos2$ + K i sin 2 6 cos 2 $ + K% cos 2 6}/2 = 

A + £cos26 + C*sin 2 6cos2$ = a + bY 2fi (&) + cY 2 , 2c (S, $) (35) 

for Ki = Hi/ 'fa where Hi play the role of Jacobi-Hooke coefficients, and 

A = K 3 /A + {K x + K 2 }/8 , B = K 3 /A - {K x + K 2 }/8 ,C={K X - K 2 }/4 , (36) 

the Y's are spherical harmonics (c and s subscripts thereon standing for cosine and sine <&-parts) and the precise form of 
the constants a, b, c is not required for this paper. This potential has as a 'very special' case B = C = 0, for which the 
potential is constant, and the 'special case' C = for which the dynamics is separable (which is sketched in Fig 3). In 
terms of the Ki, the special case corresponds to K\ — K 2 , i.e. that each cluster has the same 'constitution': the same 
Jacobi-Hooke coefficient per Jacobi cluster mass, which is a kind of 'homogeneity requirement' on the 'structure formation' 
in the cosmological analogy. The very special case then corresponds to K\ = K 2 — K 3 , for which high-symmetry situation 

7 Smith 1621 pointed out this generalization but not, as far as we are aware of, its rational interpretation. 
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the various potentials can balance out to produce the constant. Additionally the B « A perturbative regime about the 
very special case signifies Ki + K 2 << K 3 so the inter-cluster spring is a lot stronger than the intra-cluster springs, which 
in some ways is analogous to scalefactor dominance over inhomogeneous dynamics in cosmology. On the other hand, the 
C « A regime corresponds to either or both of the conditions K\ + K 2 « K3, K\ ss K 2 the latter of which signifies 
high contents homogeneity. The multiplicity of forms of writing the potential above is useful to bear in mind in searching 
for mathematical analogues for the present problem in e.g. the molecular physics literature (c.f. Sec 3.6). 
If one started instead with springs between all 6 pairs of particles, one would obtain V 6 = V + V' for 



V' = 



EL 



Uv?n k = £>sin 2 6 sin 2$ + £sin 26 cos 4> + Fsin 26 sin $ = dY 2t2s (@, $) + eF 2 ,i c (9, $) + /Y 2 ,i s (9, $) (37) 



where i, j, k are a cycle so that the first form of Vg is the most general homogeneous quadratic polynomial in the n 1 , 
D = L3/2, E — L 2 /2 and F = Li/2, and the detailed form of the constants d, e and / are not needed for this paper. 
On the face of it, this is a more general problem than the preceding paragraph's, with three further nonseparable terms. 
However, there is a sense in which these three terms can be made to go away, c.f. Sec 2.12. One can imagine whichever 
of these problems' potentials as a superposition of familiar 'orbital shaped' lumps, though such a superposition will of 
course in general alter the number, size and position of peaks and valleys according to what coefficients each harmonic 
contribution has. Contrast with the triangleland model is also interesting at this point - there Y and just two of the 
first-order spherical harmonics arose. 



The equations of motion for this potential are 



e* 



sin 6 cos 6 $* 2 = sin 26 {2B - C cos 2$ - D sin 2$} + 2 cos 26 {E cos $ + F sin $} , 



{sin 2 6 }* = 2 sin 2 6 {C sin 2$ - D cos 2$} + sin 26 {E sin $ - F cos $} , 
one of which can be replaced by the 'energy' first integral 



(38) 
(39) 



{6* 



2 6 <I>* 2 }/2 + A + 5 cos 26 + C sin 2 6 cos 2$ + D sin 2 6 sin 2$ + E sin 26 cos $ + F sin 26 sin $ = E . (40) 



Then if C = D = E = F = 0, one has a special potential, so the $ Euler-Lagrange equation gives another first integral 



( 28 ) and the subsequent quadrature for the shape of the orbit is 



$ - $ = ±T> J dO ^/ sin 6 ^2{ E - A - B cos 26}sin 2 6 - V 2 . 



(41) 





A>B>0 



-A < B < 



Figure 3: We sketch V over the sphere for the mechanically significant cases a) A > B > and b) — A < B < 0. The first has barriers at the 
poles and a well around the equator, while the second has wells at the pole and a barrier around the equator. Each is a surface of revolution 
of the curve provided, the first case being a peanut or ellipsoid-like surfaces and the second case being a concave or convex wheel. In each 
case, considering V c ff for T> 7^ adds a spike at each pole (this now means that for T> ^ both islands cannot simultaneously collapse to their 
generally distinct centre of mass points). Finally, note that our potential is axisymmetric and refiectible about its equator so its symmetry 
group is Boo x Z2 (D denotes dihaedral). If this is aligned with a DD axis of the physical interpretation, the overall problem retains a D4 X Z2 
symmetry group, of order 16. 

One simple consideration here is small and large regimes for the special case. More precisely, these are near-North 
Pole and near-South Pole regimes in 6 but become large and small regimes in terms of 1Z = tan^. For this (including 
changing to the tilded banal representation), and using 'shifted energy' E' = E — A — B 

W = E - V = 4E'/{1 + ft 2 } 2 + 32Bft 2 /{l + K 2 } 4 . (42) 
Then the near-North Pole regime (1Z << 1) maps to the problem with flat polar kinetic term and 

W = 4E' + 8{4S - E'}"K 2 (43) 



up to 0(7?. 4 ). This has the mathematics of a 2-d isotropic harmonic oscillator, 

W = £ - uj 2 K 2 /2 , (44) 

provided that the 'classical frequency' (for us with units of //time) u < (else it would be a constant potential problem 
or an upside-down harmonic oscillator problem), alongside £ > to stand a chance of then meeting classical energy 
requirements. Writing £ and ui 2 out by comparing the previous two equations, these inequalities signify that 2E > K 3 and 
2E > K3 + 2{Ks — Ki}, the latter being more stringent if K% > K\ ('stronger inter-cluster binding') and less stringent if 
K3 < K\ ('weaker inter-cluster binding'). One can also deduce from the first of these and A3 > (spring) that E > 0. 
Next, note that the near-South Pole regime (1Z << 1) maps to the problem with flat polar kinetic term and 

W = 4E'/ft 4 + 8{AB - E'} /II 6 (45) 



up to 0(1 /H 8 ) . Moreover, U — 1 /1Z maps the large case's (42 1 to the small case's ( 45 1 , so this is also an isotropic harmonic 



oscillator - in (U, <I>) coordinates and with the same £ and uj as above. One of us had previously observed a 'large-small' 
duality of this sort in triangleland [9]. It halves the required solving to understand 0«O and it. Another lesson learnt 
from the triangleland study is that we know that study of second approximations is considerably more profitable than 
that of first approximations, so we pass straight to them. Note that, for our subsequent QM study, we want the isotropic 
harmonic oscillator rather than cases corresponding to other values of the parameters £ and —cu 2 (e.g. the upside-down 
isotropic harmonic oscillator). 

2.9 Classical solutions for V = 

= V = sin 2 ©<I>* so either sin© = and one is stuck on a pole or $ is constant. In terms of the rf, this translates 
to n y — kn x for k constant, so motion is restricted to lying on a diameter. In the case of V = (w hich corresponds to 
constant potential case with T> = 0), one likewise obtains n y = kn x = ln z for I also constant, but (16 1 holds too, so all n l 



take fixed values and motion is restricted to a point. Being purely 1-d motions or 0-d motions, this subsection's solutions' 
simpleness renders them of limited interest. 1-d motions include 1) going up and down the 1-axis, corresponding to cluster 
{34} always being collapsed while cluster {12} varies in size including going through zero size at the origin and two triple 
collisions in which each of particles 1 and 2 coincide with the collapsed cluster. 2) The {12} <-> {34} of this going up and 
down the 2-axis. 3) Going up and down an n y — ±n x line, corresponding to the clusters always being of the same size 
(contents homogeneity) but that size varying from zero ({12,34} DD collision, i.e. {12} collapsing to a point and also {34} 
collapsing to a point) to maximal [in which the two clusters are superposed into the {13,24} or {14,23} DD collisions]. 

2.10 Classical solution in the very special case 

For D ^ 0, the very special case is solved by the geodesies on the shape space sphere, 

cos($ - $ ) = Kcot O (46) 



for k = V I -\/2{E — ^4} — 2? 2 , constant. Then in terms of the n % (or RelSize variables), (17 i-iv) gives 



nn z = n x cos <&o + n y sin <3>o j (47) 

i.e. restriction to a plane through the origin, with arbitrary normal (cos<I>o, sin$o; —n). But also J2 3 - 1 { nl } 2 = 1' so 
we are restricted to the intersection of the sphere and the arbitrary plane through its centre, which is clearly another 
well-known way of describing the great circles as circles within M 3 . 

The disc in the equatorial plane is particularly useful for considering the mechanics of the problem with clusters {12} 
and {34} picked out by our choice of Jacobi H-coordintes. Eliminating n z projects an ellipse onto this disc, 

k 2 = {k 2 + cos 2 $ }n x 2 + 2 cos $0 sin $ n x n y + {k 2 + sin 2 ^}^ 2 , (48) 

centred on the origin with its principal axes in general not aligned with the coordinates. E.g. for <f>o = 0, the ellipse is 

{RelSize(l,2)/{l + K ~ 2 }- 1 / 2 } 2 + RelSize(3,4) 2 = 1 , (49) 

which has major axis in the RelSize(3,4) — n y direction and minor axis in the RelSize(l,2) = n x direction, while the 



value of RelSize (12, 34) — n z around the actual curve can then be read off (47) to be n z — n x /n. With reference to the 
first subfigure in Figure 2a), as T> — > 00, the dynamical trajectory is the equator, corresponding to maximally-merged 
configurations including four DD collisions. For T> small, the motion approximately goes up and down a meridian, e.g. 
forming a basic unit of a narrow cycle from the polar DD to slightly around the T on the Greenwich meridian (reflections 
of) which is repeated various times to form the whole trajectory. [The actual limiting on-axis motion T> = is excluded 
from this subsection's working but already considered in the preceding one.] 

Other $0 straightforwardly correspond to rotated ellipses. However the mechanical meaning of these differs. E.g. about 
7r/2 clusters {12} and {34} are interchanged, while about ±7r/4 also has distinct sharp physical significance. Throughout, 
note the periodicity of the motion (already clear in the spherical model as the great circles are closed curves). The 
tessellation lines are great circles, projecting to the disc rim, the axes, the lines at 7r/4 to the axes and ellipses with 
principal directions aligned with the preceding. 
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2.11 Approximate classical solutions in the special case 

At the level of the sphere and using the second approximation, we can transcribe the solution from [9 to be, with T> for 
J and with our £ and uj 2 in place of that paper's Qo and Q2, and defining f — 2£/V 2 , / 2 = uj 2 /V 2 and g = \J /q 2 — f-2, 



V/o + gcos(2{$-$ }) = l/n = cot"- 



(50) 



In terms of TZ these are ellipses centred on the origin (including the bounding case of circles but excluding the other 
bounded case of pairs of straight lines (Fig 4) . 



circles 





ellipses 



hyperbolae strai S ht lines 



GO 



Figure 4: a) For the small approximation in (TZ, $) coordinates, the desired parameter space is the indicated wedge populated by ellipses; 
on the bounding parabola, we get circles, coinciding with the edge of the n x , n y disc for /o = 1 = /2 and becoming smaller in either direction, 
b) For the large approximation in (W, $) coordinates, the parameter space is likewise; if one converts to (7?., "3?) coordinates, however, the 
wedge is then populated by ellipse-like curves and peanut-like curves. 

N.B. we now have a third inequality on £ and uj: A£ 2 > T> 2 uj 2 , that replaces £ > as it is more stringent. Thus in terms 
of E and the Ki we get our allowed wedge of parameter space to be 2{2E — K3} 2 /T> 2 > {2E — K3} + 2{K\ — K3} > 0. 
Saturation of this corresponds to circular trajectories. For such circles to exist, the discriminant gives the condition 
{2?/4} 2 > K3 — Ki, so that the relative dilational quantity is bounded from below by the amount by which the inter- 
cluster spring dominates. 



Then by (17 v-vi) and (18 1, (50 1 becomes 

y fo + g{{n x 2 - % 2 } cos 2<J> + 2n x n y sin 2$ }/{n x 2 + n y 2 } = y/ {\ + n z }/{l - n z } , 
so solving for n z and applying the on-sphere condition for <!>n = 0, say, gives 

a 2 {/o + .9-lK 2 + {/o 
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IK! 

{fo + 9+ l}^ 2 + {fo~9+ IK 2 



(51) 



(52) 



Then one can write down a curve in terms of two independent variables such as RelSize(l,2) and RelSize(3,4): cither 
RelSize(l,2) + RelSize(3,4) = (so both are because they are positive and so both clusters have collapsed) or 

{{fo +9 + l}RelSize(l,2) 2 + {/„ - g + l}RclSize(3,4) 2 } 2 = 4{{/ + . 9 }RclSize(l,2) 2 + {/ - «?}RelSize(3, 4) 2 } . (53) 

The large regime then has 



v / /o + ffcos(2{$-$ }) = 1/W = K = tan 



2 ' 



(54) 



which is, for the cases of interest, an ellipse-like or peanut-like curve (see 9 a and Fig 4b). Applying this paper's in- 
terpretation in terms of n l or RelSize variables, the same answer as for the small regime arises again. This conclusion 
just reflects the potential imposed having an antipodal symmetry, which physically translates to shapes and their mirror 
images behaving in the same fashion. 

Finally, note that this approximate problem has a part-hidden SO(3) symmetry (such are well-known for harmonic 
oscillators). In the present context, however, its objects take the form 



Hi = um x n y +T>-sT>2/un x n v , H 2 = T> = D 2 n x /n y - T>xn y /n x , H 3 = 



2 }/2 + D! 2 /2um x 2 -D 2 2 /2ojn y 2 . (55) 



Thus its unhidden part is TL2 — 2? which has relative dilational momentum significance, while its remaining hidden parts 
are mixed shape and dilational objects. [This is used as an example in the discussion in Appendix A.] 



2.12 Discussion of 'more general combinations of springs' 

The Li terms (or, equivalently, D, E and F terms) can be dropped in the sense that one can pass to normal coordinates 
for which the symmetric matrix of Jacobi-Hooke coefficients has been diagonalized. Unlike in triangleland, however, this 
does not send one to the special case - the C-term survives and so requires addressing separately (e.g. perturbatively) . 
The elimination thus of D, E and F terms is also subject to the mechanical interpretation of the normal-coordinate 
problem being more difficult algebraically than for D = E = F = 0, so that it is conceivable that one might prefer to 
retain this simpler interpretation and treat D, E and F perturbatively. 
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3 4-stop metroland at the quantum level 



Kincmatical quantization |66] for this problem |19j involves three objects u l whose squares add up to 1 (which in the 
present case we identify with the Cartesian unit vectors n l ) and three SO(3) objects (which in the present case are the 
T>i). We then consider the time-independent Schrodinger equation [19 

{sinei^^sine^^J^+lsine}" 2 *,^ = {A + £ cos 26 + Csin 2 9 cos 2$}* , (56) 

where A = 2{A — E}/h 2 , B = 2B/H 2 , C = 2C/H 2 are dimensionless constants. Note that the above equation is separable 
for = C i.e. = C i.e. K\ — Ki\ most of our work is for this case. 



3.1 Explanation of the choice of operator ordering 

We choose an ordering that is coordinatization invariant on configuration space [67], i.e. a member of the family D 2 — 
£Ric(M) [c.f. ([8] |9|] where D 2 and Ric(M) are the Laplacian and the Ricci scalar corresponding to the kinetic metric 
M on configuration space. Moreover, following from the appropriateness of relational actions for whole-universe physics, 
observing that these have banal conformal invariance as a simple and natural feature and then asking for this to hold at 
the quantum level in the whole-universe context (i.e. in quantum cosmology or toy models thereof), among the preceding 
family of orderings we are uniquely led to the conformal ordering, for which £ = {k — 2}/4{fc — 1}, where k is the 
configuration space dimension. (See |68j for more on this motivation for conformal ordering, previous motivation on 
different premises for it being in e.g. |49l I69j .) Moreover, presently we are in configuration space dimension 2, for which 
£ = so our operator ordering choice is, in this case the same as the Laplacian ordering (itself advocated in e.g. [70] , 
while [71 also considered 2-d configuration spaces so the Laplacian conformal ordering coincidence also applies). 

3.2 Solution in very special case 



The C = case of Eq ( 56 1 separates to simple harmonic motion and the equation 

{sine}" 1 {sine* : e},e-{sine}" 2 m 2 * = ^* + 6cos26* , (57) 

If B = as well - our very special problem -, then from Sec 2.4 this has similar mathematics to ordinary QM's cen- 
tral potential problem, in which the quantum Hamiltonian H, total angular momentum Lxotai = Sq=i L q 2 anc ^ m &g- 
netic/axial/projected angular momentum L 3 form a complete set of commuting operators and as such share eigenvalues 
and eigenfunctions. In fact (also Sec 2.4) our very special problem is mathematically the same as the rigid rotor, for which 
H is Lxotai up to multiplicative and additive constants, so, effectively one has a complete set of two commuting operators, 
whose eigenvalues and eigenfunctions are the well-known spherical harmonics and, moreover also occur as a separated-out 
part of the corresponding scaled relational particle model problem. However, our 'rigid rotor' is in configuration space 
rather than in space and with total relative dilational momentum ©Total = m place of total angular momen- 

tum and projected relative dilational momentum 2?3 in place of axial angular momentum. These then have eigenvalues 
h 2 D{D + 1} and hd respectively, so we term D and d respectively the total and projected relative dilational quantum num- 
bers (which are analogous to the ordinary central force problem/rigid rotor's total and axial/ magnetic angular momentum 
quantum numbers). 

Our very special problem's time-independent Schrodinger equation separates into simple harmonic motion and the 
associated Legendre equation (in X — cos0) i.e. the spherical harmonics equations, Thus its solutions are 

* D d(6, $) cx y D d(6, $) cx Pd(cos 6)exp(±id$) (58) 

for Pp(X) the associated Legendre functions of X, D S No and d such that |d| < D. Also, D{D + 1} = —A, which, 
interpreted in terms of the original quantities of the problem, is the condition 

E' = E-if 3 /2 = /i 2 D{D + l}/2 (59) 

on the model universe's 'energy' and inter-cluster effective spring in order to have any quantum solutions (E is fixed as 
this is a whole-universe model so there is nothing external from which it could gain or lose energy) . If this is the case, 
there are then 2D + 1 solutions labelled by d (we can see the preceding sentence cuts down on a given system's solution 
space, though the more usual larger solution space still exists in the 'multiverse' sense [5U]). 

Furthermore, using a basis with sines and cosines instead of positive and negative exponentials, 

* Da > 1 ) cxAAK) . (60) 

Here, the D-label runs over the orbital types (s for D = 0, p for D — 1, d for D = 2 ...) and M is the 'naming polynomial' 
i.e. 1 for s, n x for p n ^, n x n y for d nxrly etc. (Note that the name 'z 2 ' in d nz 2 is indeed shorthand for z 2 - 1/3; shorthand 
begins to proliferate if one goes beyond the d-orbitals; the polynomials arising in our working are also subject to being 
'nonunique' under ^2^ = i{n 1 } 2 = 1.) That the wavefunctions are their own naming polynomials is via Sec 2.4's analogy 2 
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mirroring how the orbitals in space historically got their Cartesian names, and also is akin to representations |72j of the 
spherical harmonics in terms of homogeneous polynomials. Another form for the solution is[^] 

f D | d |(n l ) oc P^(n z )T d (n x /^n x 2 + n y 2 ^ = P rj (RelSize(12,34))7 d f RdSize(l,2)/^/ 1 - RclSize(12,34) 2 ^ . (61) 

However, via Sec 2.3's tessellation trick, we can interpret the wavefunctions in terms of the metroland mechanics on the 
sphere itself, on which they take the particularly familiar 'orbital' form. 

For D, d = 0, (s-orbital), note that the axis is arbitrary so it is evident from using 2 different principal axes that 
the probability distribution function on-axis is not to be trusted in spherical coordinates about that axis. We conclude 
that the ground state does not have bias toward any particular configurations. For D, d = 1, (p Uz orbital), equatorial 
configurations are improbable, meaning that mergers of the {12} and {34} clusters [including the non {12,34} DD's] are 
disfavoured, while polar configurations are probable, meaning that the {12} and {34} clusters being small and well-apart 
is favoured. For D, d = 1, ±1, in the p Hy orbital case, the n y = RelSize(3,4) axis part of the equator is probable, so 
mergers of a small {12} and a large {34} are favoured. p ri:c is the {12} <-> {34} of this. For D, d = 2, (d nr 2 orbital), 
both equatorial and polar configurations are probable, so that the {12} and {34} clusters are either merged or small and 
well-apart. For D, d = 2, ±1, in the Pn v n z case, equatorial configurations are improbable, so mergers of {12} and {34} 
are improbable, and also the {12} cluster is small; DD's are disfavoured. rf„ x „ z is the {12} <-> {34} of this. For D, d 
— 2, ±2, in the d nx 2_ ny 2 case, equatorial configurations, i.e. mergers of {12} and {34}, are probable, especially those 
with one but not both of the clusters are large (i.e. configurations along one of the RelSize (12) or RelSize(34) axes: 
contents inhomogeneity) , including the {13,24}, {14,23}, {23,14} and {24,13} DD's. In the d„ x „ case, again equatorial 
configurations are probable, but now with |RelSize(12)| ss | RelSize (3, 4) | i.e. contents homogeneity, including the DD's 
where the two clusters are on top of each other [{13,24} and {14,23}]. 

3.3 Overlap integrals: shapes and spreads of shape operators 

We are interested furthermore in computing overlap integrals (Didi |Operator|D2d2) for three applications 1) expectation 
and spread of shape operators (below). 2) Time-independent perturbation theory about the very special solution in Sec 
3.5. 3) Time-dependent perturbation theory on space of shapes with respect to a time provided by the scale in the shape- 
scale split ERPM models in semiclassical formulation also makes use of these. This parallels Halliwell-Hawking's work 
[S"2~] and embodies one of our program's eventual goals, so we prefer giving details of computing the overlaps to giving 
details of 2). 2) and 3) have the merit of extending to far more general potential terms than the harmonic oscillator- 
like terms discussed in the present working, while 2) survives as a subproblem in the corresponding time-independent 
non-semiclassically approximated shape-scale split ERPM. 

The idea to use 1) such can be traced back to how expectations and spreads of powers of r are used in the study 
of atoms (see e.g. [73] [75] for elementary use in the study of hydrogen, or [74, for use in approximate studies of larger 
atoms). Doing this amounts to acknowledging that 'modal' quantities (peaks and valleys), as read off from plots or by 
the calculus, are only part of the picture: such as the mean (nlm | r | nlm), (nlm | r 2 | nlm) and the spread A n i m (r) = 
■\/(nlm|r 2 |nlm) — (nlm|r|nlm) 2 are also useful. E.g. for hydrogen, one obtains from the angular factors of the 
integrals trivially cancelling and orthogonality and recurrence relation properties of Laguerre polynomials in Appendix B 
for the radial factors that 

(nlm|r |nlm) = {3n 2 -1{1+ l}}a/2 and A nlm r = ^{n 2 {n 2 + 2} - {1{1 + l}} 2 }a/2 , (62) 

where a is the Bohr radius. One can then infer from this that a minimal typical size is 3a/2 and that the radius and 
its spread both become large for large quantum numbers, c.f. how the modal estimate of minimal typical size is a 
itself; the slight disagreement between these is some indication of the limited accuracy to which either estimate should 
be trusted. Also, we identify the above as expectations of scale operators, and thereby next ask whether they have pure 
shape counterparts in the standard atomic context. 

The answer is yes. Up to normalization, they are the 3-Y integrals |75j (for Y spherical harmonics, the radial parts of 
the integration now trivially cancelling) , and the general case of this has been evaluated in terms of Wigner 3j symbols 
|75j . Furthermore, many of the integrals for the present paper's specific cases of interest are written out case-by case in 
[7B"] (this applies to expectations of the RelSize 's as well as B, C, D, E, F perturbation terms' constituent overlaps). 
Shape operators for hydrogen are also considered in [77] (briefly) and |78j . 

Moreover, the context in which shape operators occur in molecular physics is wider than just the above. 
E.g. 1) expectations of cos/3 for (3 a relative angle from inner products between physically meaningful vectors e.g. between 
the 2 electron-nucleus relative position vectors in Helium, in the characterization of molecules' bonds or in nuclear 
spin-spin coupling (p 443 of [77]). 



E.g. 2) one also gets expectations of Y2q(0) [c.f. form 4 of (35 1] in spin-spin and hyperfine interactions (p 437-441) of [77] 
(as a shape factor occurring alongside a 1/r 3 scale factor. 



8 This is found by shifting from arctan to arccos and then using one of the standard definitions of Tchebychev polynomials, T d (£) = 
cos(d arccos(£)) Despite being the product of two generally nonpolynomial factors, the two conspire to produce polynomials in each case. We 
then introduce the symbol 7d(£) to mean T d (£) for cosine solutions and y/l — T d (^) 2 for sine solutions. 
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E.g. 3) In the study of the H£ molecular ion, one uses fixed nuclear separation as a scale setter and then one has not 
only 1 relative angle but also 2 ratios forming spheroidal coordinates with respect to which this problem separates, and 
expectations of all these things then make good sense. 

We contemplate 'mini- and midi'superspace counterparts of such shape operators in Sec 4.5. 

As regards good shape operators for 4-stop metroland, the kinematical quantization carries guarantees that the three 
n l are promoted to good quantum operators. These can be interpreted as RelSize(l,2), RelSize(3,4) and RelSize(12,34) 
as per Sec 2.3. It is also useful to note at this stage that n z is not only physically RelSize(12,34) but also mathematically 
the Legendre variable. 

Then (D d | RelSrzefl,2) | D d) = (D d | RelSize(3,4) | D d) = and (D d | 6 | Dd) w (D d | RelSl^e(l2,34) | D d) = as an 
obvious result of orientational symmetry. The useful information starts with the spreads, 



A Dd (RelSi Z e(l,2)) = f { ^^ } QM , A Dd (RelSM3,4)) = J^^^Q^) , (63) 



A Dd (9) « A Dd (RelSi Z e(12,34)) = \j ^^-Ij^+S} 1 ' (64) 

for (52(d) = 1/2 for the d cosine solution, 3/2 for the d sine solution, and 1 otherwise, and Qi(d) the sin <-> cos of this. 
One can then readily check that (n x 2 + n y 2 + n 2 ) = 1, as it should be. 

One case of interest is the ground state. Therein, the spreads in each are l/VS. Another case of interest is the large 
quantum number limit. Ao d (0) ss Arj d (RelSize(12,34)) which, for the maximal d (|d| = D), is equal to l/\/2D + 3 which 
goes as l/\/2D — > for D large. The hydrogen counterpart of this result is AuO ~ l/v2i — > 0, i.e. restriction to the 
Kepler-Coulomb plane (e.g. [77] outlines this, while [78] considers it in more detail). Back to our problem, this result 
therefore signifies recovery of the equatorial classical geodesic as the limit of an ever-thinner belt in the limit of large 
maximal projectional relative dilational quantum number |d| = D ('the rim of the disc' of Sec 2.10, traversed in either 
direction according to the sign of d). In fact, as for the constant potential we can put the axis wherever we please, this 

leads to recovery of any of the classical geodesies. Also, for d = 0, Ado© ~ Ad o(RelSize(12,34)) — ► 1/V2 for D large. 
This means that the s, p riz , d nr 2 ... sequence of orbitals does not get much narrower as D increases, so that for these 
states we only get limited peaking about clusters {12} and {34} both being small and well apart, a situation which we 
will revisit in the next subsection due to its centrality to the assumptions made in, and applications of, this paper. The 
RelSize(l,2) and RelSize(3,4) operators' spreads tend to finite constant values for large D no matter what value d takes. 

What of <&? Now, clearly, by factorization and cancellation of the O-integrals, the d = states obey the uniform 
distribution over to 2ir, with mean it and variance 7r 2 /3 (corresponding to axisymmetry) . Furthermore, (D d | <i> | D d) is 
also 7r and cosine and sine states have 

r2 /Q J_ 1 /OA 2 anrl Ar. *{tf>\ — */tt2/o — 1 /OA 2 



A Dd ($) = V^/S + l/2d" and A Dd ($) = ^/tt 2 /3 - l/2d 2 , (65) 

which indicate some resemblance to the uniform distribution arising for large d (mean and variance do not see the 
multimodality, but at least, by inspection along the lines of the preceding subsection, it is regular multimodality for d 
maximal - equatorial flowers of 2D petals - by inspection of the shapes of the standard maximal s, p, d, f, g ... orbitals. 

3.4 Solution in special case — large and small regimes 

Passing to stereographic coordinates, banal-conformal transforming to the flat representation and applying the small 
approximation, our Schrodinger equation becomes 

- {n 2 /2}{Rr 1 {m >n }^ i + n^ 2 ^^} = £- lu 2 K 2 /2 , (66) 

which is in direct correspondence with the 2-d quantum isotropic harmonic oscillator (see e.g. [73] [TS] [EH] under 1Z < — ► r 
(radial coordinate), 1 < — > particle mass, and with our oj as classical frequency (x/). Thereby, 

£ = Uhuj for n = 1 + 2n + |d| (67) 

for n a node-counting quantum number running over No and d a 'projected' dilational quantum number as in the preceding 
subsection but now running over Z. [The 'shifted energy' in its usual units, E' — E — A — B, itself goes as 

E' = {n 2 ^ 2 /2}{l + y/l - B{A/\lh} 2 } , (68) 

so for Uh/uj << 1 (small quantum numbers as used below), E' I w UhH, for f2 = 2y/—B.] The solutions are then (to 
suitable approximation) 

f Nd (e,$) cx 6 |d| {l + |d|e 2 /12}exp(-cje 2 /8S)Lj, d| (we 2 /4S)exp(±id$) (69) 
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for L a (£) the associated Laguerre polynomials in £ (see Appendix B). [The <!>-factor of this is rewriteable as before in 
terms of the n % or RelSize(12,34) and RelSize(l,2), while the 0-factor is now a somewhat more complicated function of 
RelSize(12,34)]. 



The large regime gives the same eigenvalue condition ( 67 1 , and ( 69 1 again for wavefunctions except that one now uses 
the supplementary angle S = tt — G in place of O. Next, see Fig 5 for the form and interpretation of the wavefunctions. 



DD 




N= 1, d = 



N — 0, Idl — 2 



Figure 5: Probability density functions for this subsection's problem for ui/h large, 400, say, plotted using Maple |81| . All d = states 
are axisymmetric about the {12,34} clustering's DD axis, i.e. all relative sizes for cluster {12} and for cluster {34} are equally favoured. The 
ground state is peaked around the {12,34} DD collision. It it the surface of revolution of the given curve. The N = 0, |d| = 1 solutions are a 
degenerate pair. Each takes the form of a pair of inclined lobes - the cosine solution's oriented about the n x = RelSize(l,2) = D collision 
and the sine solution's about the n y = RelSize(3,4) = D collision. These next three solutions form a degenerate triplet. The N = 1, d = 
solution is a slender bulge around the {12,34} DD collision, then a gap and then a second bulge in the form of a cone, representing 'a band 
very close to this DD collision and a band somewhat close to it being probable, while all other configurations are improbable. The N = 0, 
|d| = 2 solutions are tulips of four petals, the cosine one separately favouring the lunes of near {12} D collisions and near {34} D collisions 
(i.e. contents inhomogeneity) , while the sine one disfavours these and favours instead the lunes at 7r/4 to the preceding, which correspond to 
contents homogeneity of the {12} and {34} clusters. 

The large case's approximate solution is just the reflection of the preceding about the equatorial plane with the same interpretation except 
that {34} is now to the left of {12}. 

In the small regime, the RelSize(l,2) and RelSize(3,4) operators still have zero expectation as each sign for these remains 
equally probable. For D, d substantially smaller than u/h so powers of the latter dominate powers of the former (and 
uj/h was considered to be large, so this works for the kind of quantum numbers in this subsection's specific calculations), 
the following mean and spread results for shape operators are derived using orthogonality of, and a recurrence relation 
for, Laguerre polynomials, as provided in Appendix B. 



(Nd | RclSize(12,34) | Nd) = 1 - 2Uh/u 



(70) 



A N dRelSize(12,34) is zero to first two orders, beyond which the approximations used begin to break down, but it would 
appear to have leading term proportional to h/uj. These results signify that the potential has trapped what was much 
more uniform in Sec 3.3 into a narrow area around the {12,34} DD collision. Furthermore, 



A Nd ^« y/2UHQa(d)/u 



(71) 



for a = 1,2 gives A N d(RelSize(l,2)) for a — 1 and A N d(RelSize(3,4)) for a — 2. So one can obtain strong concentration 
around the poles by suitable choice of springs, amounting to a tall thick equatorial barrier and polar wells. The ground 
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state has the tightest spread in RelSize(l,2) and RelSize(3,4): \j2hjio. This has some parallels with how the Bohr radius 
is an indicator of atomic size, including the hydrogen-isotropic harmonic oscillator correspondence |79j . 



3.5 Perturbations about the very special solution 

We begin by recasting our Schrodinger equation in Legendre variables n z — cos0, 

{{1 -rc z 2 }# n J,„ 8 + {l-n z 2 }- 1 ^^ = {A-B + n z 2 {2B~Ccos2^}+Ccos2^>}^ . (72) 

One can then study this using time-independent perturbation theory (see e.g. |75j for derivation of the formulae for this 
up to second order). Applying perturbation theory here means considering 1) C small, which is high contents homogeneity 
at the level of each cluster's (Hooke cooefficient)/(reduced mass) in the sense that K\ — K 2 is small compared to ti 2 . 2) 
B small, in the sense that h 2 is large compared to {Ki + K 2 }/2 — K 3 , which collapses to Ki — K 3 small in the case of 
C — 0, meaning that there is little difference between the inter-cluster spring and the intra-clustcr springs. 



Perturbative study of ( 72 ) is amenable to exact calculations though involving various of trigonometric and stan- 
dard/tabulated associated Legendre function integrals, or, alternatively, the aforementioned 3-Y integrals. Furthermore, 
this continues to be the case if one includes a non-diagonal/non-normal basis' D, E and F terms. 

For the S-perturbation, as both it and the unperturbed Hamiltonian commute with T>, the eigenvalue problem can 
be solved separately in each subspace V<j of a given eigenvalue d of V, and in each such subspace the spectrum of 
the unperturbed Hamiltonian is nondegenerate, so that nondegenerate perturbation theory is applicable (this argument 
parallels e.g. p 697 of [73] ). This gives (with the unperturbed problem's A playing the role usually ascribed to the energy 

and Hi the perturbative term) = (D d | H x | D d) at first order and A^\ =-J2t>> ,a' ?n ,d I ( D ' d 'l H i I D d ) l 2 /{-V - A u } 
at second order [75]. Then e.g. [19] double use of a standard recurrence relation [82] gives a Ad = 0, AD = 0, ±2 'selection 
rule'. Moreover, the terms that survive this take the following forms. 

(Dd\B{2n z 2 - l}|Dd> = £{l-4d 2 }/{2D-l}{2D + 3} , (73) 

which is closely related to the expectation of n z — RclSize(12,34) already computed in Sec 3.3. Two new overlaps that 
are more general than expectations are 

(D+2d | B ^- 1}|Dd) ^/ «^^;»g:;f-^ , 

and then, swapping D for D - 2, also, 



(D - 2 d | B{2n 2 - 1} | D d) = -** ^ ~ ~ ~ ^ . (75) 

V 1 x z s 1 l 2D - iy {2D + 1}{2D - 3} V ' 

Using these then gives the perturbed 'energies': 

E Dd = A + fi 2 D{D + l}/2 + B{\ - 4d 2 }/{2D - 1}{2D + 3}+ 

4£ 2 {{2D + 5}{2D + 3} 3 {D 2 - d 2 }{{^ - l} 2 - d 2 } - {2D - 1} 3 {2D - 3}{{D + 2} 2 - d 2 }{{D + I} 2 - d 2 }}} 3 

fi 2 {2D + 5}{2D + 3} 3 {2D + l}{2D-l} 3 {2D-3} 1 >' 

(76) 

Note that d positive and negative are treated the same, so there is only a partial uplifting of degeneracy. Changes to the 
wavefunction due to the perturbations for the sign of B corresponding to Sec 3.4 and to second order in B are that we 
get slight bulges at the poles for the ground state (a bit of d nz 2 mixed in). 

The C perturbation can likewise be studied based on half-way stage overlaps that can be directly transcribed by our 
angular momentum to dilational momentum analogy from those computed in e.g. |76j . E.g. the surviving terms are found 
[19] by a second standard recurrence relation [82] to obey the selection rule Ad = ±2, AD = 0. Some noteworthy features 
of the study of the C term are that degenerate perturbation theory is now required, there is no first order contribution as 
Ad = ±2 only, and now d and -d do get shifted differently corresponding to this perturbation not preserving the axis of 
symmetry. In nondiagonal/nonnormal form, the further D term has the same selection rule to the C term's while the E 
and F terms share the selection rule Ad = ±1, AD = 0, ±2. The above 'noteworthy features' apply to these also. 

3.6 Molecular physics analogies 



Analogy A) (57 1 occurs in mathematical physics (e.g. from the separation of the wave equation in prolate spherical 
coordinates [831 1541 [551 152"] ) and has multiple applications in molecular physics studies of which parallel some of the 
studies in the present paper. Examples of this in molecular physics are as follows. 



Analogy A.l) (57) recast in terms of the Legendre variable is 

{{l-n z 2 }*, n J, nz -{l-n z 2 }- 1 mH = {A-B+2Bn z 2 }y , (77) 
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which is the easier of the two spheroidal equations that arise in the study of the HJ molecular ion [86j |87l [88] . This and 
the next two analogies are for B < 0, although the aforementioned mathematical physics literature covers B > 0. 



Analogy A. 2) The potential Vb{l — cos 28} [c.f. form 3 of ( 35 1] occurs in modelling the rotation of a linear molecule in 
a crystal [5U1 [55]. Here, the analogy is ( 21|20|22" l where the axis and rotor in question are provided by the linear 



molecule itself, 'energy' <-> energy up to a constant, 

Ki/2 <-» 2Vq up to the same constant difference as in the energy analogy (78) 

B «-> -2V (79) 



Analogy A. 3) The potential — Q||E 2 cos 2 6 | [c.f.form 2 of ( 35 1] for ot\\ the polarizability along the axis occurs in the study 
[9"2l |9"3"] of e.g. the CO2 molecule in a background electric field E (the study of polarizability is the theory underlying 
Raman spectroscopy). Here the analogy is, rather, 

B < — ► -a|,E 2 /2 . (80) 

Analogy A. 4) Examples 2) of Sec 3.3 is another substantially developed area in the molecular physics literature. 
Analogy B) is with the ammonia molecule NH3, in the following rougher but qualitatively valuable sense. NH3 has two 
potential wells separated by a barrier and then is capable of tunnelling between the two at the quantum level (like an 
umbrella inverting in the wind). Our model for B < is similar to this, albeit in spherical polar coordinates: we have 2 
polar wells with an equatorial barrier in between. 

This analogy then gives us some idea about how the separate solutions for the two wells compose. For NH 3 , one can start 
with separate solutions for each well and additional degeneracies ensue (due to the wells being identical and being able 
to distribute some fixed energies between these in diverse ways). However the wavefunctions tend to perturb each other 
toward breaking these degeneracies, forming symmetric and antisymmetric wavefunctions over the two wells |77L 187] . 

3.7 Applications of the analogies and developing an overall picture of our model 

Firstly, the B < locally stable small or large regimes are the kind of regimes that are termed 'rotator-like' in analogy 
A.2)'s literature; both of the SO(3) quantum numbers (for us, dilational quantum numbers) hold good in this regime. 

Secondly, for B < one can use Analogy B to form a simple picture of putting the small and large approximations 
together. As d remains a good quantum number for the unapproximated problem, one expects to need the North Pole 
approximation's d and the South Pole approximation's d to match and the subsequent perturbations exacted by these 
two approximations upon each other not to affect d. Also, prior to any recombination, one has degeneracies as follows 
(call the near-North Pole's node-counting quantum number n and the near-South Pole's n'). There is the one ground 
state n = n' = d = 0, then the degenerate pair n = n' = 0,d = ±1, and then the degenerate quadruplet n = 1,n' = 
or n = 0, n' = 1 for each of d = ±1. Now if n and n' match, expectation of RelSize(12,34) goes to again though 
the wavefunction's distribution is bimodal about both poles. If they do not match, RelSize(12,34) retains some nonzero 
expectation due to the peaking near the two poles being different in detail. The flip here, as in NH3, is an inversion, i.e. 
it reverses the orientation, sending 1,2,3,4 to 4,3,2,1. 

Thirdly, analogy A. 3) is well-known for its Raman-type ±2 and not ±1 selection rule, which parallels our results of 
Sec 3.5. Analogy A. 3) has furthermore been studied perturbatively for what for us is the small B < regime. This allows 
us to e.g. check the half-way house results ([73] [74] [75]) against p 271-273 of [9"3"] . 

Fourthly, further resources from analogy A.l)'s references [HH [55] include analysis of this equation's poles in the 
complex plane and how it admits a solution in the form of an infinite series in associated Legendre functions in the 
vicinity of ±1 and in Bcsscl functions in the vicinity of 00, as well as how to piece together these different representations. 
It is then appropriate to compare results from the expansion in associated Legendre functions against our perturbative 



regime (this particular working holds regardless of the sign of B). Thus we find the lowest four cases of ( 76 ) to agree with 
p 1502-4 of [55], which additionally provides the corresponding wavefunctions which we use to first order in B in Sec 4.4 
to evaluate the naive Schrodingcr interpretation probabilities for these states' model universes being large. 

Fifthly, one cannot really put together our near-polar calculations and our perturbative calculations, because the "B 
small perturbative condition" goes a long way toward u> being small and then only a bit of the wavefunction is near the 
pole. Our near-polar calculations should be compared, rather, with the asymptotics for B large. Analogy A.2)'s literature 
covers this for what for us is the B large negative (>> E — A — B = E') regime, giving, via the analogy, 

E' ~ n/kj largc + 0(1 /^largc) and (81) 

* cx exp( Wlargo cosO/ft){{tanf } 2N {secf } 2 {l d l +1 > + 0(l/w largo )} (82) 

as the relevant asymptotic solutions, for wi argc = 2\J—B. Now, from ( |43|44 1 the small-0 approximate solution (69 1 's 
uj = 8^—B = 4wi argc , so exp(wi arge cos6/fi) in (81 1 « const x exp({w/4/i}{— 6 2 /2}), which is indeed in agreement with 



the leading and dominant factor of (69). For our model, this regime signifies that A'3 >> K\,Kq, i.e. that the inter- 



cluster spring is much stronger than each of the intra-cluster springs. This is termed a 'harmonic oscillator-like regime' 



- comparing (81) and the standard result for the 2-d isotropic harmonic oscillator makes it clear why. d alone is a good 
dilational quantum number in this regime. 

Sixthly, the spheroidal equation has led to many hundreds of pages of tabulations |84j and further numerical work e.g. 
in jB2][M], though the most recent of this states that this study is still open in some aspects. 

Seventhly one can furthermore envisage extending analogy A. 2) to have the further parallel with our model that a 
rotationally-dislocated molecule in a cubic crystal will have preferred directions in space of approximately the same form 
as ours are in configuration space. We do not know if such a study has been done. 

Finally we comment that A. 3) has been extended [55] to include what for us are C, D, E and F terms. For, what 
one has more generally is a symmetric polarization tensor a such that fi p — a pa E a . Then for the CO2 model in a 
diagonal basis a z = a\\ giving the combination -a± sin 2 # — a\ \ cos 2 # [a slight improvement of analogy A. 3) by inclusion 
of the smaller a± = a x = a y ], and this readily rearranges to the special case of the third form of (351. But for more 
general groups than just oxygen atoms at each end of the axis (while still remaining in a diagon al b asis) a x a y , 
giving — a x sin 2 6*cos 2 — a y sin 2 #cos 2 </> — a z cos 2 # which is the general case of the second form of (35). Moreover, in 
non-diagonal bases, the off-diagonal elements form extra terms directly analogous to those in (37). Thus there is an 
extended analogy between our problem and the study of polarization, with 4-stop metroland's Jacobi-Hooke coefficients 
forming a configuration space-indexed analogue of the spatial-indexed polarizability tensor. 



4 Conclusion 

Relational particle models (RPM's) benefit from notions of locality and structure that are absent in minisuperspace and 
are free of many of the technical difficulties of midisuperspace models. This makes them suitable for testing some of the 
conceptual aspects of quantum cosmology, and of quantum general relativity (such as the Problem of Time) . In particular, 
in this paper we study the RPM of 4 particles in 1-d - 4-stop metroland - in the case without scale, both classically and 
quantum-mcchanically. We concentrate on the clustering into two particular binary clusters [of particles {12} and particles 
{34}] both by using coordinates that follow this case and imposing a potential term that restricts the physics to being 
near such a configuration. This is toward a qualitative conceptual model of the quantum cosmological seeding of structure 
formation in a semiclassical regime (paralleling the Halliwcll- Hawking 52J approach, which is somewhat narrower as a 
Problem of Time strategy but has further conceptual and computational applications outside of the Problem of Time 
context too), and of records theory [13] [55J EE] IH [55] [57] . The counterpart of the current paper's model with scale (which 
is harder and in which the current paper's work occurs as a subworking under the shape-scale split) will be required 
for some aspects of such a study (in particular, for a semiclassical treatment with a greater number of parallels to that 
of GR). This is further work in progress [TO] HH [12] [28], though Sec 4.2-4 give a brief account of generalizations of the 
current paper's model and how these meet additional quantum cosmological and Problem of Time criteria. 

This paper's model has an S 2 configuration space and then the mathematics which follows has analogies with the 
standard axisymmetric sphere and central force problems of ordinary mechanics. In particular, where a conserved angular 
momentum occurs in these analogue problems, a conserved relative dilational momentum occurs in our model. [These 
both have SO (3) mathematics, but each has a different physical nature, the two being embraced by our notion of rational 
momentum which generalizes angular momentum to ratios that do not happen to physically be angles.] We then interpret 
some of 4-stop metroland's classical and quantum solutions in cases with harmonic oscillator-like potentials. The solutions 
in spherical variables give fairly standard mathematics such as that of the rigid rotor and of the 2-d isotropic harmonic 
oscillator in some of the simpler cases, albeit now these require subsequent interesting and unusual interpretation in 
terms of the 4-stop metroland problem's mechanical variables. We deduce this at the level of mass-weighted coordinates 
by tessellating the shape space sphere by the mechanical interpretation appropriate to 4-stop metroland, which we find 
to possess the symmetry group of the cube. Further tools we introduce, paralleling basic treatises on the atom, are 
expectations and spreads of shape operators, to which we can also attribute cosmological analogies. Our shape operators 
are RelSize(12,34): the relative size difference between universe and its {12}, {34} cluster contents, RelSize(l,2): the size 
of the {12} cluster relative to the size of the whole model universe, and its {34} cluster counterpart. The polar angle $ 
itself is an inhomogcncity ratio of the contents of the universe themselves (i.e. of the two clusters relative to each other). 

We obtain expectations and spreads for such operators e.g. in ground state and in large quantum number limits. We 
consider the very special constant potential case that resides within the harmonic oscillator-like potential models as a 
particularly structurally homogeneous balance of springs, as well as more general cases treated perturbatively for small 
differences in spring constitution, asymptotically for large such differences, and in near-polar approximations. We further 
benefit from recognizing that the special case with the two clusters of the same spring constitution but the inter-cluster 
spring is weaker gives an elsewise well-known spheroidal equation, alongside various molecular physics analogies: with , 
NH 3 , rotation of molecules in crystals and molecular polarizability (at least the last of which extends to cases with more 
general combinations of springs between the four particles). This permits us to tap into substantial mathematical physics 
results and generally control our particular problem, and is a further example of useful bridges between RPM quantum 
cosmology models and the physics of molecules (triangleland RPM with harmonic oscillator-like potentials having already 
been found to share mathematics with the Stark effect for a linear rigid rotor [19]). 
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4.1 Comments on extension to N > 4 metrolands 



The Jacobi H- and K-coordinates of Fig 1 generalize to cover an increasing variety of 'part H-shaped, part K-shaped' 
clusterings may be useful in this respect), which are of additional value as less trival models of structure formation and 
of records theory. For full reduction for scalefree arbitrary- N-stop metroland and the subsequent Euler-Lagrange equations 
in the arbitrary-potential case, see [18(19]. Moreover, we now comment that the number and nature of conserved quantities 
that each of these possesses is tied to the usual SO(N - 1) representation theory, and their physical interpretations extend 
the present paper's discovery of dilational quantities. Tessellations by physical interpretation are now harder as they both 
have more pieces and also are more difficult to visualize due to being higher-dimensional. However, our relative size and 
contents inhomogeneity shape operators do straightforwardly extend to N-stop metroland. 

Within each N-stop metroland, one can envisage a tower of special, very special, ... (very) N ~ 2 special problems. The 
most special of these in each case has a constant potential and thus gives ultraspherical geodesies classically and the 
ultrasphcrical rigid rotor quantum- mechanically (solved by ultraspherical harmonics |19j). while the next most special of 
these in each case has (N - l)-d isotropic harmonic oscillator mathematics in its near-polar regime (solved by a power 
times a Gaussian times an associated Laguerre polynomial). Establishing a perturbative regime about each most special 
problem would then appear to be possible e.g. [19] by recurrence relations of the Gegenbauer polynomials [82j|97j. One 
technical difference is that, if one does use conformal operator ordering, then one can no longer use the configuration 
space being 2-d to evoke collapse to Laplacian ordering like in Sec 3.1. However, hyperspheres are of constant curvature 
and so of constant Ricci scalar curvature, so £Ric(S fc ) is just a constant, ^k{k — 1} (our 'E has no nonconstant prefactors 
banal conformal representation' having the unit sphere as its configuration space). So, even in this case, the sole difference 
between Laplace and conformal ordering (or any other member of the D 2 — £Ric(M) family of operators) is in what is to 
be interpreted to be the zero of the energy. We also note that for N = 5 the analogy with the Halliwell-Hawking scheme is 
somewhat tighter, as both involve perturbative expansions in § 3 ultrasphcrical harmonics. Finally, the next most special 
equation unapproximated can also be mapped to the spheroidal equation, so that the fairly standard mathematical physics 
of that equation continues to be of aid in N-stop metroland. 

4.2 Comments on Extension to metrolands with scale 

The configuration spaces for these are cones over the corresponding shape spaces [TT]. The present paper's advances in 
the physical understanding of conserved quantities in RPM's have further applications here. Our introduction of shape 
quantities to be promoted to operators also continues to be relevant here through there being a shape-scale split, so that 
evaluating shape operators here collapses back to pure shape workings such as the present paper's. Expectation and 
spread of size (in close parallel with atomic physics) will also now be pertinent. Metrolands with scale will now have 
solid rather than surface analogues of the present paper's 'tessellation by physical interpretation' technique. The way in 
which similarity RPM arises as a subproblem from the shape-scale split of the scaled theory means that the most special 
harmonic oscillator case and perturbations thereabout survives as a piece of the analysis upon introduction of scale, now 
partnered by isotropic harmonic oscillators in the size quantity. 

This setting with scale is more appropriate as regards both toy-modelling of cosmology in general and in particular 
to using a semiclassical approach both in the Problem of Time context and in the Halliwell-Hawking context. One of us 
makes a first sketch at this in the smallest case in [10 , with other cases to follow in [TT1 IT21 128] . 

4.3 Comments on extension to 2-d 

The price to pay in introducing scale is that one no longer has nontrivial constraints in spatially l-d models with which to 
model some effects due to GR's momentum constraint. On the long term, one can get around this by passing to spatial 
dimension > 1. A first such model is triangleland: the RPM of 3 particles in the plane. This has a S 2 shape space like 
the present paper's 4-stop metroland model does. This gives a number of useful insights. E.g. parts of the present paper 
parallel • Even more significantly, because the 4-stop metroland interpretation of the sphere turns out to be more 
straightforward, the present paper allows for an improved understanding of the less straightforward triangleland case [20] . 
Scaled triangleland is harder; so far we have just provided some classical study for this (also to be augmented by the 
present paper's techniques at the classical level in [11] toward finally providing a quantum study of it in |14j). 

Moreover, triangleland lacks the present paper's nice feature of splitting into two nontrivial subsystems (of 2 particles 
each), which is a useful nontriviality from the structure formation and records theory perspectives. Studying scaled 
quadrilateralland (RPM of 4 particles in the plane) [98 would incorporate this feature too. Thus this model would 
possess a number of midisuperspace's features with the benefit of being technically simpler. This makes it particularly 
suitable for the simultaneous investigation of records theory and the semiclassical approach (which may support each 
other, and histories theory, to form a more robust combined approach to the Problem of Time and to quantum cosmology 
[9"9l 1551 [H)]). Quadrilateralland does have a further technical complexity - its shape space is CP 2 , which unavoidably 
involves complex-projective mathematics (triangleland has CP 1 but this is well-known to also be § 2 ). 

Further features for consideration in RPM models involve [13] A) oriented shapes - real projective spaces MP N ~ 2 in 
place of § N ~ 2 as shape spaces or CP N ~ 2 /Z 2 in place of CP N_2 as shape spaces and the corresponding cones in models 
with scale and/or B) (partial) particle indistinguishibility by which only pieces of whichever of the preceding spaces would 



pass to being the configuration spaces. The present paper's treatment of physical interpretation by multiple coordinate 
charts and tessellations by physical interpretation are doubtlessly ideas of further value in the study of these models with 
their wide range of configuration space geometries. 

4.4 Further details of Problem of Time applications 

Our wavefunctions, eigenvalues and operators are useful in the following Problem of Time investigations. 

1) The computation of nai've Schrodinger interpretation |53j probabilities of the universe having some particular property. 
Example 1) consider quantifying P(universe is large), in the sense that the two clusters under study are but specks in the 
firmament, by P(e-close to the {12,34} double-double collision), which means, at the level of the configurations themselves, 
that the magnitude of -\/RelSize(l, 2) 2 + RelSize(3, 4) 2 /Relsize(12,34) lies between 1 and 1 - e 2 /2, and, in configuration 
space terms, that one is in the e-caps about each pole. Then from the latter and by the naive Schrodinger interpretation, 
this probability oc J e _ ca s j^ r | 2 dS' = f^ZoUe=o + J^L^JW©; $)| 2 sin 6 d6d<i>. So, e.g. for the very special solution's 
ground state and first excited state, one gets proportionality to e 2 + 0(e 4 ), while for the states with dilational quantum 
numbers D = 1, |d| = 1, one gets proportionality to e 4 + 0(e 6 ). 

Example 2) consider quantifying P(the two clusters nominally under study are in fact merged) by P(£-close to {12,34} 
merger) which means, at the level of the configurations themselves, that the size of Relsize(12,34) does not exceed the 
small number 5, and, in configuration space terms, that one is in the <5-belt around the equator. Thus 
P(<5-close to {12,34} merger) oc / 5 _ bc i t | ^ 1 2 dS', which, in the very special case, works out to be proportional to <5 3 + 0(<5 5 ) 
for D = 1 d = and to S + O(S) 3 for the other three lowest-lying states. 

Example 3) consider quantifying P(universe is contents-homogeneous) in the sense that the two clusters under study are 
similar to each other, by the magnitude of RclSizc(l,2)/RelSize(3,4) departing from 1 by no more than 2r\. Then, on 
configuration space, one is in the tetralune described in Fig 2, and the naive Schrodinger interpretation gives 
P(universe is ^-contents-homogeneous) oc /^_ totra i unc |^| 2 d5, which, in the very special case, comes out as proportional to 
r\ for all four of the lowest-lying states. 

Example 1 also makes sense for the small-regime special solution. One now obtains proportionality to e 2 y/ui/H to 
leading order i.e. the same '(small) 2 ' factor as in the very special problem but now with an opposing '-^/large' factor, 
amounting to the small regime's potential well (Fig 3b) concentrating the wavefunction near the poles i.e. in the region 
of the configuration space corresponding to large universes in the above-described sense. 

Finally, also repeating Example 1 for the wavefunctions with first order perturbative corrections in B included, we now 
find proportionality to e 2 {l - 8BI 2 /9h 2 } + 0(B 2 ) + 0(e 4 ) for the ground state, to e 2 {l - 8BI 2 /25h 2 } + 0(B 2 ) + 0(e 4 ) 
for D = 1, d = 0, and e 4 {l - 36I 2 B/25H 2 } + 0(B 2 ) + 0(e 6 ) for the D = 1, d = 1 states. The signs of these corrections 
conform with intuition, as (Fig 3a) B > corresponds to placing a potential barrier at the poles and a well around the 
equator, which should indeed decrease the amount of wavefunction there, i.e. making large universes less probable, and 
vice versa for B < 0. 

2) Given explicit wavefunctions such as this paper's, one can build up projectors and mixed states (including with 
environment portions traced out) and then construct conditional probabilities |43j for pairs of universe properties. 

3) As regards records theory |57j . the current paper's classical work provides some means of defining a notion of distance on 
configuration space (which is quite closely related to the measure problem in cosmology [H]), and a notion of localizability 
in space. Next |100j . one would construct notions of information (alias negentropy) both at the classical level and at the 
quantum level for the problems solved in this paper. In this respect it is worth noting that QM perturbation theory 
suffices in order to build an approximate statistical mechanics [101J. Such notions of information include e.g. Shannon's, 
von Neumann's, Tsallis's [57], as well as notions of subsystem information, mutual information and correlation (such as 
the covariance for the two clusters in the situation that the present paper centres on). 

4) One can also build up decoherence functionals for histories theory [58], and consider the Feynman- Vernon influence 
functional that Halliwcll uses [56] for the study of records within histories theory. 

5) Exact wavefunctions also serve as useful checks on whether the semiclassical approach's assumptions and approximations 
are appropriate [TH H5J . 

N.B. Operator insertions for meaningful shape operators remain useful in constructing various of these Problem of Time- 
relevant objects. All of 2) to 5) above workings being substantially longer than 1), we leave them for future occasions. 

4.5 Analogues of our shape operators in mini and midisuperspace? 

Another longer-term goal would be to export insights acquired by our program to 'mini and midi'superspace; are there 
then useful analogues of shape operators for these (anisotropy operators, inhomogeneity operators?) In surveying the 
literature, we have found, firstly, that Kuchaf and Ryan |102j consider (y 2 ) for y a reparametrization of one of the 
anisotropy degrees of freedom f3± in diagonal Bianchi IX quantum cosmology. Secondly, Ashtekar and Bojowald make 
mention of an anisotropy operator in studying loop quantum gravity 103 . Thirdly, Petryk and Schleich [104 consider 
expectation values for geometrical quantities in the Hartle-Hawking initial state in their study of conditional probabilities 
in the 3-d Ponzano-Regge minisuperspace. Fourthly, Halliwcll and Hawking [52 compute the expectation of the anisotropy 
in temperature of the microwave background; this has the additional value of being "halfway to midisuperspace" in that 
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it considers inhomogeneous perturbations about a homogeneous spacctimc. As regards inhomogeneous spacctimes, the 
Lemaitre-Tolman-Bondi solution principally concerns radial scale variables, so it is far more of an analogue to the shape- 
scale extension of the present paper. While there is not anything as yet that we know about using shape operators at the 
quantum level for the Gowdy universe, e.g. Andersson, van Elst and Uggla [ 105) use a form of shape-scale variables at the 
classical level. Thus, while shape quantities and shape operators have occasionally been used in 'mini and midi'superspace, 
a systematic treatment parallelling that in atomic physics does as yet appear to be lacking. 

Acknowledgments: we thank Claire Anderson for hospitality, and Professors Don Page, Jonathan Halliwcll and Julian 
Barbour for references and discussion. This research was partly supported by Grant Number RFP2-08-05 from The 
Foundational Questions Institute (fqxi.org). 

Appendix A This paper's generalization of angular momentum put into context 



What is habitually called angular momentum mathematics (because that is a common guise in which it appears in 
physics), is, de facto, the representation theory of SO(p). The most usual case is that of SO(3), habitually called the 
rotation group, and also interpretable as the isometry group of the 2-sphere [more generally, SO(p) is the p-dimensional 
rotation group and the isometry group of the (p - l)-sphere]. In turn, SO(3) is closely related to SU(2) (which is its double 
cover). The rational momentum viewpoint is more general than the mechanical angular momentum perspective but not 
the SO(p) one, as Fig 6 begins to explain by classifying examples. Beyond that, one can get SO(p) as a piece of an even 
larger group. [This is already the case for Fig 6's Runge-Lenz example, which can be viewed as a second SO(3) coming 
from a partly-'hidden' SO(4), but it covers further cases where the group is not necessarily 'hidden', such as SO(3, 1) in 
relativistic particle physics, or SU(3) containing 3 different directions' worth of SU(2) ladder operators]. Quantum 'SO(p) 
objects' also combine under an addition rule, whereby composites of cases in Fig 6 arise. E.g. total angular momentum 
T = L + S in atomic physics; more generally, total rational momentum T = 1Z + A, which also would include the effect 
of adding internal 'arrows' to the present paper's 1-d RPM setting]^] There are also rational momenta that are linear 
combinations of relative angular momentum and relative dilational momenta occur in triangleland [20 . Finally, there are 
SO(p) quantities that are physically phase space objects. In this sense, this paper's rational momentum came about from 
building a different notion of 1-d space from the usual into the present paper (scalcfrccncss) Finally, adding arrows to the 
present model could be nontrivial in the sense that the 1-d arrows need not just obey separate 'tensored-on' occupation 
rules since dilational momentum exists as does addition of arrow and dilational momentum quantities, so that spatially 
1-d models can have arrow-dilational momentum interactions that parallel the spin-orbital angular momentum couplings 
that occur in higher spatial dimensions. 
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Fi gure 6; Various physical realizations of SO(p) objects. 

Appendix B Some results concerning associated Laguerre Polynomials 

The bounded solutions of the associated Laguerre equation 

+ {a + i- Ov,z + Py = o 

are the associated Laguerre polynomials L^(£). These obey [32] the orthogonality relation 

/>oo 

/ rexp(-£)L^)L^)d£ = unless = (3 1 
Jo 



and the recurrence relation 



£L%(0 = {2p + a+ 1}L|(0 - {f3 + 1}L£ +1 (0 - {/? + a}L^(0 



(83) 
(84) 

(85) 



As the word 'spin' itself has rotational and hence angular momentum connotations, we generalize 'spin alias internal angular momentum', 
S, to 'arrow alias internal rational momentum', A. 
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